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It has been proposed that antiferromagnetic Fe adatom spins on semiconductor Cu-N surfaces 
can be used to store information [S. Loth et al, Science 335, 196 (2012)]. Here, we investigate 
spin dynamics of such antiferromagnetic systems through Monte Carlo simulations. We find out 
the temperature and size laws of switching rates of Neel states and show that the Neel states can 
become stable enough for the information storage when the number of spins reaches to one or two 
dozens of the Fe spins. We also explore promising methods for manipulating the Neel states. These 
could help realize information storage with such antiferromagnetic spin systems. 

PACS numbers: 75.75.-c, 75.78.-n, 75.10.-b, 75.90.+W 


I. INTRODUCTION 


II. SPIN MODEL AND SIMULATION METHOD 


The atomic-size antiferromagnets (AFM) on semicon¬ 
ductor surfaces have been realized thanks to technology 
advance, and have attracted a lot of attention for their 
potential in magnetic storage technology[Il-l5|. Due to in¬ 
sensitivity to magnetic fields, these nanomagnets are con¬ 
sidered to be able to meet the needs of the vast growth 
of storage density in magnetic media 0. On the other 
hand, it is also far more difficult to effectively manipu¬ 
late the AFM stagger magnetization (or Neel state) than 
ferromagnetic magnetization in ferromagnetic materials. 
Modern scanning tunneling microscopy (STM) technol¬ 
ogy can be used to manipulate magnetic states of such 
nanomagnets through applying an electric field or imect- 
ing an electron current on one of adatom spins [5l-l9|. It 
has been reported that AFM chains assembled by plac¬ 
ing Fe adatoms on a CU 2 N overlayer on Cu(lOO) can be 
switched between two quasistable Neel states when the 
polarized electrons are made to flow across one atom of 
the chainQ. There have been theoretical studies on mag- 
netic transitions induced by tunneling electrons |lfll - ll4l |. 
but temperature effects need to be considered. A reli¬ 
able study on the temperature-dependent dynamics of 
Neel states is in need. 

Here, we use the dynamic Monte-Carlo method to 
study the switching rates of a series of such antiferro¬ 
magnetic Fe spin chains. We investigate effects of vari¬ 
ous temperatures and spin-polarized currents in terms of 
experimental parameters. Arrhenius-like behaviors are 
observed until the temperature is too low to hurdle the 
thermal activation barrier of the spin reversal. Our re¬ 
sults show that such spin chains can be stable enough for 
practical usage when one or two dozen Fe adatoms are 
used. We also explore effective methods for external ma¬ 
nipulation of the Neel states. More detailed results will 
be presented in the following. 
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For low temperatures (T < 10 K), the Fe spin bi-chains 
on the CuN/Cu(100) surface can be effectively described 
with an anisotropic Heisenberg antiferromagnetic model 
of^ = 20,i, 


^0 = E 

where 


( 1 ) 


= -giisB ■ kj - - [Sl^f] ( 2 ) 

describes the f-th Fe spin in the j-th chain (j = 1, 2). 

Here, Sij = {Sf is the spin vector operator 
for the single Fe spin, and as shown in Fig. 1, J is the 
nearest antiferromagnetic exchange constant in one of the 
chains, and J' the ferromagnetic exchange constant be¬ 
tween the two chains. Letting J' = 0 and confining j to 
1 in (1), we can obtain a single chain system. The pa¬ 
rameter g is the Lande g factor (here g = 2 \s used), and 
fiB the Bohr magneton. The parameters D and E are 
used to characterize the single-ion magnetic anisotropy, 
and B = is the applied magnetic field. In 

addition to (1), an effective magnetic field Bup needs to 

be applied on the hrst spin Sip to split the two Neel 
states and some spin polarized electron currents must 

be injected through S'lp to control and detect the Neel 
states, as have been done in the previous experimental 
and theoretical papers [3, [3|. 

We use the theoretical methods by Gauyacq et alfl^ 
to treat the effects of injected electrons on the spins 
in the chain and bi-chain systems. When we treat a 
spin at the (i,j) site to calculate the thermal activa¬ 
tion barrier, we use mean-field approximation for its 

nearest spins, Si+ij = Si+ij. This results in effec¬ 
tive Hamiltonian Hq = where Hij = —Sij ■ 

F),, - D(%)2 - E[{SN)^ - with F)., defined as 
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FIG. 1. (a) A schematic of Fe adatom spins on a CU 2 N over¬ 
layer on Cu(lOO), with electronic current injected through 
the STM tip. (b) The inter-spin exchange constant along the 
chain is J. (c) The inter-spin exchange constant between the 
two chains is J'. 


(b) 

*e 
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JSi+ij + J'Sij+i + gfisB. Here, S^+ij and S^j+i are 
already classical quantities, with the transverse compo¬ 
nents being equivalent to zero. Their starting values are 
set in terms of the Neel states (more detail will be given 
in the following), and they are updated at each step of 
the Monte Carlo simulation. For each electron tunneling 
from the STM tip through the first Fe atom and into the 
substrate, there is a chance to change the states of the 
whole Fe spin chains. 

For our DMC simulation, the MC steps are defined 
by the time points: = At • n, where n takes non¬ 

negative integers in sequence. The At is set to 1.6 x IQ-^ 
s or determined by At = le/J, where / is the current 
intensity from the STM tips. At is chosen to satisfy the 
requirement that there is only one electron within a MC 
step. At the beginning, we set all of the spins at one 
for Neel state Eq, and the Ath spin should be set as 
(—1)-^ • S. For the nth step, the MC simulations are 
performed in the following way. 

For each Fe spin, we take Hi j as a 5 states quantum 
model, and there is a thermal-activated energy barrier 
Ae between the starting and ending states. Using the 
Arrhenius law[I^ [ 1 ^ , we obtain the thermal-activated 
spin reversal probability, Pt = 1 — exp(rAt), for each 
MC steo pTHl^ . where r = roexp(-^^) is the spin re¬ 
versal rate and rp is the character attempt frequency. 
When an injected electron is tunneling through the first 
Fe spin, the collision channel states Sc are defined to 
simulate the electron-Fe-spin coupling, and the electron- 
activated spin reversal probability is defined as Pe BUS- 




FIG. 2. Temperature dependence of the switching rates with¬ 
out magnetic field for 1 x A linear chains with A = 4, 6, 8, 
and 10 (a); and 2 x A bi-chains with A = 4, 6, 8, and 10 (b). 


[^ . Therefore, the total probability for the first Fe spin 
is equivalent to P/ = 1 — (1 — Pe)(l — Pt). Actually, at 
extra-low temperature an additional collective spin chan¬ 
nel becomes available for the system to transit from one 
Neel state to the other. This transition probability, calcu¬ 
lated through exact diagonalization of the Hamiltonian, 
is added to the above probability to completely describe 
the spin dynamics. 

For the nth MC step, each spin in the system has 
chance to reverse. With n increasing, at last, the sys¬ 
tem finally transits to the other Neel state Ei at a spe¬ 
cial n. This n value is denoted by A^.. We define 
Ti = (Nr) ■ At as the average transition time from Eq to 
El- The average transition time from Ei to Eq, T2, is cal¬ 
culated in the same way. The switching rate is defined as 
R = 2/{T 1 +T 2 ). To be consistent with experiment 0, |^, 
we take J=-0.737 meV, J'=0.03 meV, P>=1.87 meV, 
P=0.061 meV, — 0.115 tesla, and ro=2xl0®s-^ in 
the following. Our main results are calculated by aver¬ 
aging over 100,000 runs of DMC simulations. 


III. INTRINSIC LIFE TIMES OF NEEL STATES 

We simulate the spin dynamics for the different Fe 
chains with magnetic fields. It is confirmed that the mag¬ 
netic field has little effect on switching rates. Simulated 
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switching rates for 1 x Fe-spin chains are presented 
in Fig. mja). It is clear that the switching rates show 
strong temperature dependence as expected. Above 6 
K, the switching rates of the 1 x iV Fe chains mainly 
follow the Arrhenius law. With the addition of two Fe 
spins, the energy barrier increases with 4.15 meV. The 
switching rates of Neel states in long Fe chains are tiny in 
comparison with those in short ones. When temperature 
decreases away from the Arrhenius regime, the switching 
rates become independent of temperature. Such behav¬ 
ior can be attributed to the quantum tunneling because 
the thermal activity is frozen out at those ultra-low tem¬ 
peratures. Simulated switching rates of 2 x iV chains are 
shown in Fig. EJb). In these cases, the Arrhenius law 
will work well until the temperature is lower than about 
3 K, and the quantum tunneling rates become so small 
that the Neel states can be considered to be quite stable 
below about 2 K. It should be pointed that the simu¬ 
lated temperature-rate curves of single chains for 1x6 
and 1x8 are in good agreement with experimental curves 
and those of bi-chains for 2 x 4 and 2x6 are also con¬ 
sistent with experimental results Q, which show that our 
methods and results are both reliable. 

It is very interesting that the temperature-dependent 
switching rates can be fitted with the simple function: 


R = Ro exp( 


-AE 

ksT 


) + Rt, 


(3) 




FIG. 3. The energy barrier (AE) and the residual switching 
rates (Rt) (a), and the life times (r) (b) of the 1 x A Fe linear 
chains, with N taking even numbers. 


where Rq, AE, and Rt are htting parameters. The in¬ 
trinsic life time of the Neel state can be described with 
T = 1/R. For T = 0, we get r = tq = 1/Rt, the low 
temperature limit of the life time. For the 1 x N chains, 
since i?o is nearly constant, Rq ^ 2.0 x 10“® per second, 
we present AE, Rt, and r (for T=0, 8, lOK) as func¬ 
tions of chain length in Fig. 3. Simple fitting leads to 
the following important results. 

AA = 2.07-A-4.37, tq = 2.97 x 10“^° x 27.9^, (4) 

where AE is in meV and rg in second (s). The tunneling 
rate Rt will reduce by 778-fold when two Fe adatoms 
are addeT which is in good agreement with experimen¬ 
tal valueQ- With all the three parameters known, we 
can calculate life times (r = 1/R) of Neel states for any 
given temperature and chain length. For the ultra-low 
temperature limit (OK), we have 0.12 s, 1.8 minute, 1.0 
day, and 2.1 year for the 1 x N chains with N=6, 8, 
10, and 12, respectively. For the 2 x N chains, we have 
Rq = 2.0 ^ 2.2 X 10“® per second, and the N dependence 
relations of AE and tq (Rt) are similar to Equ. (4). For 
a special case of 2 x A {N = 8) chains, we already have 
To ^ 28.2 day. Although increasing temperature will de¬ 
crease the life time, very long life time r of a few years 
can be achieved as long as N is large enough, making 
1 X 14 and 2 x 8 or longer chains. 


IV. MANIPULATION OF NEEL STATES 

A current injected through the STM tip on the first 
adatom spin can change the Neel states in the Fe spin 
chains Il3|. To characterize a Neel state, we define a 
Neel weight Wn of the 1 x A chain to be 
with i running from 1 to A. It is A and —A for the 
Neel states Eq and Ei, respectively. We present the Neel 
weight Wn of 1 X 8 spin chain under different current 
intensities and temperatures in Fig. |4Ka). It can be seen 
that the Neel weight changes exponentially with time. 
The temperature and the spin-polarized current play the 
key role in determining the switching rates and target 
Neel weights, but the initial status of Neel weights, rjo, 
do not affect the target Neel weights {Woo)- Generally, 
we can fit the time-dependent Neel weights with a simple 
function: 

Wn = Wo exp(-t/to) + Woo- (5) 

The to parameter describes the time scale, and Wq + Woo 
is equivalent to the initial Neel weight. The calculated 
three fitting parameters are compared in Table (jll be¬ 
tween different temperatures, different currents, and dif¬ 
ferent initial Neel weights. The temperature dependence 
of the target Neel weights for four currents (1, 2, 4, and 
8 pA) is presented in Fig. SKb). The current dependence 
of the target Neel weights for different temperatures (6, 
8, and 10 K) is presented in Fig. SKc). It is clear that 
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FIG. 4. Time dependence of the Neel weights {Wn) for the 
1 x 8 chain (a), temperature dependence for different cur¬ 
rent intensity (1, 2, 4, and 8 pA) (b) and current dependence 
for different temperature ( 6 , 8 , 10 K) (c) of the target Neel 
weights (Woo) in the same chain. 




Time(lO's) 

FIG. 5. Current manipulation of Neel states (parameterized 
with Wn) under T = 6 K for the 1 x chains with N—A (a) 
and A =8 (b,c). The current is switched between 8 pA and -8 
pA for (a) and (b), and between 5.8 and 6.2 pA for (c). 


TABLE I. The different fitting parameters for different cur¬ 
rents, temperatures, and initial status. 



r(K) 

Wo(s-i) 

to(s) 

IFoc(s-i) 

7=4pA, ryo=4 

6 

12.24 

59.77 

-4.24 


8 

12.86 

4.01 

-4.85 

7=8pA, r;o=4 

6 

15.98 

82.35 

-7.96 


8 

15.95 

6.57 

-7.95 

7=4pA, 770 =-4 

6 

-3.74 

59.78 

-4.25 


8 

-3.16 

4.01 

-4.85 


cyclic periodic function. Such manipulation of the 1 x A 
chains {N = 4, 8) with two currents with different time 
periods is presented in Fig. [S](a) and (b). The Neel states 
are almost fully switched for these two cases. The two 
curves look similar although the time scales are hugely 
different from each other. Furthermore, the big step in 
the I-Woo curves in Fig. 4(c) can be used to manipu¬ 
late the Neel states. When the current is periodically 
switched between 6.2 pA and 5.8 pA, the Neel weight is 
periodically switched between 5.1 and 7.9, as shown in 
see Fig. HKc). 


the temperature dependence is weak, but the target Neel 
weight increases substantially when the current intensity 
increases. It is interesting that there exists a big Woo step 
at ^6pA. It is because one more spin channel becomes 
available at the equivalent voltage for inelastic tunneling 
of the electrons (24 . 

The sensitive current dependence of the target Neel 
weights can be used to manipulate the Neel states 
through suitable spin-polarized currents and tempera¬ 
tures. When a periodic spin-polarized current between 
-8 and 8 pA is applied, the Neel weights can assume a 


V. DISCUSSION AND CONCLUSION 

It was pointed out by Gauyacq et a/[l^ that there are 
three mechanisms for electron-induced switching between 
the Neel states of Fe antiferromagnetic chains. At finite 
temperature, thermal effects will go into action and in¬ 
crease exponentially with temperature. They can be at¬ 
tributed to a barrier-hurdling thermal activation mecha¬ 
nism described by an Arrhenius law. 

In summary, we have investigated spin dynamics of 
the antiferromagnetic spin systems through Monte Carlo 
simulations. We have found out the temperature and size 
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laws of switching rates of Neel states for such spin sys¬ 
tems. Furthermore, we have shown that the Neel states 
can be made stable enough for the information storage 
if one or two dozens of such Fe spins are used. We also 
have demonstrated promising methods for manipulating 
the Neel states. We believe that these can be useful to 
realize antiferromagnetic information storage. 
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